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In this paper we continue the work begun in 2002 on the identification of the analyt- 
ical expressions of Feynman integrals which require the evaluation of multiple elliptic 
integrals. We rewrite and simplify the analytical expression of the 3-loop self-mass inte- 
gral with three equal masses and on-shell external momentum. We collect and analyze a 
number of results on double and triple elliptic integrals. By using very high-precision nu- 
merical fits, for the first time we are able to identify a very compact analytical expression 
for the 4-loop on-shell self-mass integral with 4 equal masses, that is one of the master 
integrals of the 4-loop electron g-2. Moreover, we fit the analytical expressions of some 
integrals which appear in lattice perturbation theory, and in particular the 4-dimensional 
generalized Watson integral. 
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1. Introduction 

Analytical expressions of many Feynman diagrams contain polylogarithmic func- 
tions of various kinds (Nielsen polylogarithms. Harmonic polylogarithms^, Har- 
monic sums^ 3, multiple zeta values, Euler sums®, etc. . . ). But there exist Feynman 
integrals which cannot be described only in terms of polylogarithms. 

At two-loop level, the discontinuity of the off-shell massive "sunrise" diagram 
with different masses is expressed by elliptic functions. At three (or more) loop level 
the situation worsens. These diagrams contain nested multiple elliptic integrals; 
the current mathematical knowledge of such integrals is scarce or missing. At this 
preliminary stage, "experimental mathematics" is the best tool. In other words, 
high-precision numerical values of the integrals are fitted^ with various candidate 
analytical expressions until an agreement is found. The equalities are then checked 
up to hundreds or thousands of digits. In this way the right analytical expression 
is identified beyond any reasonable doubt. This may be the starting point of the 
subseque nt sear ch for a rigorous proof of the result, task which may take months of 
hard worljniIB] 

In Ref. [5] we work out a very high-precision value of the 3-loop scalar master 
integral of the "sunrise" diagram S3 of Fig.l; then we fit that value with products 
of elliptic integrals, checking the equality with a precision of thousand of digits. 

In this paper we continue the work on the 3-loop integral, and we simplify the 
3-loop analytic result by using some identities between elliptic integrals. Next we 
review and develop the approach used for fitting and identifying the 3-loop integral 
and we apply it to the 4-loop scalar master integral of the "sunrise" diagram Si of 
Fig.l. The analytical calculation of this master integral is also of physical interest, 
because it is the simplest non-trivial master integral of the 4-loop electron g-2 
(of which a high-precision numerical calculation is under way^l). Very likely, the 
analytical constants which appear in the expression of S4 should also appear in the 
analytical expression of 4-loop g-2. We calculate an high-precision value of £4 and 
we are able to fit this value with an expression containing two new elliptic constants, 
checking the equality with a precision of thousand of digits. 

We apply this procedure also to the values of some 4-dimensional lattice inte- 
grals, and we identify their analytical expressions; surprisingly, they contain the 
same elliptic constants of the 4-loop integral. 

The plan of the paper is the following: In section [2] we simplify the results of 
Ref. [5]by using identities between elliptic integrals. In section|3]we study the 4-loop 
"sunrise" integral. We collect a number of results on a "simplified" version of the 
integrals involved. Then we use these results as a guide for identifying the candidate 
analytical expressions suitable for fitting the 4-loop results. In section |4] we show 
the analytical results found for the 4-loop integrals. In section [5] we fit the values of 
some 4-dimensional lattice integrals with the same analytical constants discovered 
in the 4-loop integrals. In section [6] we give our conclusions. 
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Fig. 1. Three-loop and four-loop self-mass diagrams. 



2. Three-loop single-scale self-mass integral 
2.1. The restttts 0/ #e/. H 

In Ref. [5] we considered the Feynman diagram S 3 {p 2 , wif , rn|, m|, m|, D) with equal 
masses mj = 1, and on-shell external momentum (see Fig.l) 

D r ,A \JD 



Sa(— lj 1, 1) 1) 1, D) = [ , 2 . , 2 -T7 £- - 
J Oi + l)(ffi + l)(8a + l 



[d^gi] [d"g 2 ] [d%] 



)((P- ?i - 92 - 93) 2 + 1) 



where 



[«*°?] = 



d B g 



7r D / 2 r 3 



-1. (i) 



(2) 



By using an hyperspherical representation for the integral, we found that the value 
of S3 could be expressed as a sum of various double elliptic integrals, the simplest 
being 



A, 



dl 



dm 



R(l,-l,-T)J R(m, I, -T)R(m, -1,-1) 

2.641 379 476 074 689 431 349. . . , (3) 



R(x, y, z) = \/ x 1 + y 2 + z 2 — 2xy — 2xz — 2yz 



(4) 



We were not able to calculate A3 in analytical form. Therefore we evaluated it at 
very high precision and we tried to fit the numerical value with various kinds of 
analytical expressions. In Ref. [5] we found that 



A3 = K(w-)K(w+) 



w± 



z± 



z±-l ' 



z± = -(2-V3) 4 (4±v / 15) 2 , (5) 
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where K is the first of the two elliptic integrals 

, N f 1 dt , N f 1 dt VI - mt 2 
K(m) = / _____ , E{m) = / v _ _ . (6) 



jo vT - Vi - mt 2 Jo Vi - ? 

We were able to fit the values of S3 in 2 and 4 dimensions Q 

47T 

S 3 (-l,l,l,lA,D=2) = ^=K(w-)K(w+) , (7) 

Via 



, 22 _ 2 577 _j 4tt /35 



S3 -1, 1, 1, 1, 1, D = 4 - 2e) = 2e- J + — + — e" 1 + _ 

3 36 Vl5 



+ ^( W _)X( W+ ) - I (i?(l - «,_)£(! - «,+) + 5E(w-)E(w+))\ + ^ + 0(e). 



(8) 

Eq.©, Eq.© and Eq.© were checked up to 30000, 40000 and 1200 digits, respec- 
tively. 

2.2. New relations between elliptic integrals 

Now we note that the arguments w± of the elliptic integrals are singular values. In 
the context of elliptic integrals k r is a called singular value if 

where r is an integer or a rational number. The arguments w± of elliptic integrals 
are singular values for r = 15 and r = 5/3 respectively, 



10_ = fci 5 , U>+ = fc 5/3 



that is 



X(l- fcl5)=v - j K(l k 5{3 ) = J, (iQ) 



(11) 



X(fc 15 ) ' if(fc 5/3 ) V 3 

K(k 5/3 ) _ Vl5- V3 



AT(A;ib) 2 

The values of elliptic integrals of second kind of Eq.© are obtained following Ref.fTOl 
E ( k r) = A r-ln ^ + K ( k r) (l ~ ^] > (12) 



4^¥K(k r ) K " V V^ 
£(l-fc r ) = ^— +J ftT(fc r )a r . (13) 

(14) 



a The two-dimensional value J7} was calculated together with Jg), but not published in Ref. 131 
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By using the values 



VT5-V5-1 y/TE-VE+1 

"15 = , «5/3 = g j (15) 



where 

and expressing K(k 5 / 3 ) by using Eq. fTTj) . we are able to rewrite Eq.([7]) and Eq.© 
in the very compact form 

1,1, 1,^ = 2) = ^, (18) 

5 3 (-l, 1, 1, 1, 1, D = 4-2e) = 2e -a + 22 e - 2+ 5H e - 1 + 6191 _ 14V^_ V5 
3V ' ' ' ' ' 7 3 36 216 P 900 w 

(19) 

Eqs. (fTU)l - p5ll were also shown by David Broadhurst in his beautiful talk given 
in Bielefeld^. After the talk, our unpublished results (|17[) - (flH|) were shown to 
David Broadhurst. 

2.3. The path to Eq. (fj) 

For sake of completeness we recall here some unpublished observations that sug- 
gested us the form of Eq.©. First we calculated analytically the simplest double 
elliptic integral: 

dx f 1 dy 



o VI - x 2 Jo y 2 \/l - x 2 y 2 



Mi 



3.437 592 909. . . , (20) 



which factorizes in the square of the elliptic integral K. 

Then we observed that in the diagram S$(— 1, 1, 1, 1, 1) the value of the square of 
the external momentum p 2 is — 1, which is different from the threshold (p 2 = —16) 
or the pseudothresholds (p 2 — —4,0). We expected that the analytical expression 
is much simpler for the on-threshold diagram than for the off-threshold diagram. 
So we considered the above graph with one mass changed: S^(— 1, 1, 1, 1, 2). Now 
the value of p 2 = —1 is on a pseudothrcshold (which are at p 2 = — 1, —9, —25). The 
integral analogous to A% is 

A' 3 = — -- -/ — - = 1.474 585 992... . (21) 

J R(l,-l,-l)J R(m,l,-4)R(m,-l,-l) y 1 
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Subsequently, as we expected the form of ^3 to be more complicated than A'3, we 
tried also products of K with different arguments, and we found Eq.©. 

3. Four-loop single-scale self-mass integral 

Now we consider the single-scale 4- loop self- mass diagram S4 (p 2 , m\ , , m|, m 2 , D) 
of Fig.l, in the case of equal masses rrij — 1 and on shell external momentum 
p 2 = —1. This diagram has 3 master integrals, the simplest being 



S 4 (-l,l,l,l,l,l,£) = 

[d D q x ] [d D q 2 ] [d D q 3 ] [d D q 4 ] 



(q 2 + l)(q 2 + l)(q 2 + \){q 2 + l)((p- q x - q 2 - q 3 - q 4 ) 2 + 1) 



(23) 



■1?, 

As already said in the introduction, this is also one of the several master integrals 
which appear in the calculation of 4-loop g-2. 

The first observation is that, fortunately, the value of p 2 = — 1 in S4 is already on 
a pseudothreshold (which arep 2 = —1,-9,-25); therefore to simplify the analytical 
structure of the integral we have not to use the mass change of the previous section 

El 



3.1. High-precision numerical values 

We need an high-precision numerical value of this integral in order to obtain a 
meaningful fit. This is worked out by using the difference equation method presented 
in Refs. [13l HU Summarizing, one raises to n one denominator of Eq. ([23|) 

[d D qi] [d D q 2 ] [d D q 3 ] [d D q,] 

(q 2 + ir(q 2 + l)(q 2 + l)(q 2 + l)((p - Ql - q 2 - q 3 - q 4 f + 1) " 1 > 

The function X^n) satisfies the fourth-order difference equation 



PiX 4 (n + 3)+p 2 X 4 (n + 2)+p 3 X 4i (n + 1) +piX 4 (n) + p 5 X 4 (n - 1) = 

24(.D-2) 4 J 3 (l)J(n) , (25) 

where 

pi = -768n(n + l)(n + 2)(n - 2D + 5) , 

p 2 = 128n(n + l)(lln 2 + (63 - 26£>)n + 11D 2 - 57D + 76) , 

p 3 = 4n(-129n 3 + (294L> - 588)n 2 + (-148.D 2 + 592D - 567)n 

+ 8D 3 - 48L> 2 + A6D + 36) , (26) 
p 4 = 2(-60?i 4 + (294D - 468)n 3 + (-485D 2 + 1499Z? - 1164)n 2 + (308L> 3 

- 1363Z) 2 + 2015Z) - 996)n - 60L> 4 + 326L> 3 - 658L> 2 + 584L> - 192) , 
p 5 = -(n -D + l)(n -2D + 3)(2n - 3D + 4)(2n -5D + 8) , 
and J is the one-loop integral 

J(n) = I • (27) 

J (9i 2 + l)" 
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Eq. (|25p contains in the r.h.s. the integral obtained from Si by contracting one line, 
which factorizes into the product of 4 one-loop tadpoles. The solution of Eq. ([25|) 
compatible with the large-n boundary condition X^(n) oc n~ L ) I 2 can be written as 
Xi(n) = C 1 X^°(n) + C 2 X^ (n)+X NH (n). The functions Xf° and X NH 

are respectively the two solutions of the homogeneous equation compatible with 
the above large-n behaviour and one particular solution of the nonhomogeneous 
equation Eq. (|25|) . The constants C\ and C2 are obtained from the 3- loop self-mass 
integrals belonging to the diagram obtained from S 4 by deleting one line, that is S3. 
The amount of calculations needed to work out and solve the systems of difference 
equations is high, so that the calculations have been performed by means of an 
automatic tool, the program SYS described in Ref. [13] 
In two dimension one finds 

5 4 (-l, 1, 1, 1, 1, 1, D = 2) = 40.2451219019305821264798187417 . . . (28) 

and in the limit D — > 4 

5 45 4255 106147 , . 

S 4 (-l, 1, 1,1, 1, 1,D = 4 - 2e) = - - 3 - ^ - — (29) 

- 141.72215618664768694996791- 521. 14654568600250441775466c 

- 3347.9933650782886117865341e 2 - 17951. 3774774809944931097622e 3 

- 101753.8165331173182139560386e 4 + C>(e 5 ) . (30) 

All the above numerical constants were calculated with a precision of over 2400 
digits; for the sake of space we show here only the first 25. 



3.2. Triple elliptic integrals 

By using an hyperspherical representation for the integral, Si(— 1, 1, 1, 1, 1, 1, D = 2) 
and the finite part of S^— 1, 1, 1, 1, 1, 1, D = 4— 2e) contain triple elliptic integrals, 
the simplest being 

f°° dl f°° dm f°° dr 

A 4 = ' 



R(l, -1, -1) Jo R(m,l,-l)J R(m,r,-l)R{r,-l,-l) 

8.749 361 490.... (31) 

We prefer to redistribute the arguments of the R functions and consider the similar 
integral 

dl r°° dm r{Vl-V^) 2 dr 



R(l,-l,-l)J R(m,-l,-l)J R(r, I, m)R{r, -1,-1) 

i8.749 361 490. . . = iA A , (32) 
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and the companion integral 

dl f°° dm /-(Vz+Vm) 2 dr 



B 4 = 



R(l, -1,-1) J R(m, -1,-1) y(vr-VSi) a R(r, I, m)R(r, -1,-1) 

9.607 815 129..., (33) 

where r± = (^/l ± y/rnj are the two zeroes of R(r, I, m). A4 and B4 are the 4-loop 
integrals analogues of the 3- loop integral A3 , and we have to find their analytical 
expressions. 

3.3. Simplifying the problem 

First of all we consider the simplest triple elliptic integral with structure similar to 
Eq.flM]): 

A= f 1 - 1 t= t I' dz = 4.335 593 665. . . . (34) 

Jo Vl-a; 2 io \Jl—y 2 Jo vl — z 2 ^Jl — x 2 y 2 z 2 

By changing one limit of integration over z to the zero of 1 — x 2 y 2 z 2 we obtain a 
companion integral analogous to Eq. 



B = f 1 * r * r y f = 6.997 563 016. . . . 

Jo vl — x 2 Jq ^/l — y 2 J\ V 1 z 2 — ly 1 - x 2 y 2 z 2 

(35) 

We expect that the study of the simpler constants A and B can help us understand 
the analytical expressions of A4 and B4. Integrating over z 

A = C -fi= t -^=K{xV) , (36) 
Jo v 1 - x 2 Jo \]l - y 2 

1 dx f 1 dy 2 2\ 



B = 7T=f/ -^==i^(l - a; V) . (37) 
Jo Vl - x 1 Jo a/1 - y 2 

Integrating over y and x 

b= [ 1 ^ K (i^^) K (i_±VTEZ) . (39) 



vT 



or 



No analytical expression is known for 4.F3 ^ 2 2^ 2 j 1^ . Trying to understand the 
reason of that, we study the following family of integrals: 



1 / 1\ Q /1 X/9 



dt K m {t) K n (l - t) (-^1 I -^=== ) (40) 
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We consider here only the integrals which have results containing elliptic constants. 
At the value m + n = 1 there is the integral 

\t =2K>(1) , (41) 

equivalent to Eq. ([2H)l : note that it factorizes in a square of if (1/2) = T 2 (l/4)/ (4-^-). 
At the value m + n = 2 we find numerically that the six integrals 

dt — jj- = B , (42) 
*^E = 2 S , (43) 



s: 

Jo 



1 K 2 (t) 

dt * W = 4A (44) 

Vt Jo VT=t 

1 K(t)K(l-t) , s 

dt y V - = 2B , 46) 

can be expressed in terms of A and B. At the value m + n = 3 wc find a surprise: 
six integrals factorizes into the fourth power of i£"(l/2). 

" dt K*(t) = il<* (j) , (47) 



/ 

■A) 



Vt 5 V2 



d t^m=^a), (49 ) 



(50) 



jT 1 * K 2 (i)K(i-t) = ^ 4 Q 



A ^» = V fi ) , (51) 



The analytical results (|42 |) -([52 | have been fitted numerically and checked up to 200 
digits of precision. We find results factorized at level 2 and 4, but not at level 3 (odd). 
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This behaviour reminds us of the non-factorization of values of Riemann ^-function 
at odd integers, and suggest to consider the constants A and B as irreducible objects. 
We can relate A and B to multidimensional £-like quadruple series. Let us consider 
the integrals 

VT=t\ K{t) ) 

These integrals correspond to the integrals (|43|) . (j45|) . (|42|) . for m = 0, 1, 2, respec- 
tively, and their values are Iq = 2B, Ii = 2A, I 2 = B. We apply the change of 
variable 

q = exp(-nK(l - t)/K{t)) or equivalently 1 - t = {6 4 {q)/6 3 (q) f , (54) 
where 6j(q) are the Jacobi Theta Functions, then 

I m =^- m J o dq (el(q)e 3 (q)^-e 3 (q) - 0i(?)0 4 (g)jU(g)) (-\ogq) m . (55) 

Expanding in series the 9 functions, and integrating over q term-by-term, I m be- 
comes a quadruple series 

DO OO OO OO , ~.y +j(h 2 ,'2\ 

u-n^Y. E E E (56) 

i= — oo j — — oc k— — oo l — — oo / 



for m = 2 the series converges, so that 



oo oo oc oo 



?'— — oo j — — oo k— — oo I — — oo 

where the prime means that the origin i = j = k = l = must be excluded in the 
summation. 



3.4. Integrals of products of homogeneous solutions 

Coming back to the 4-loop integral, if we close the 4-loop self-mass diagram S4 by 
connecting together the two external lines, we obtain a 5-loop vacuum diagram. The 
5-loop vacuum diagram can be decomposed into two 2-loop self-mass diagrams con- 
nected together. Therefore in D — 2 dimensions its value is given by the integral of 
the square of the 2-loop self-mass diagram J cPp Sf (p 2 , 1, 1, 1). The vacuum diagram 
is expected to have the same analytical structure of S4, but with higher transcen- 
dentality. S^Cu) satisfies a second order nonhomogeneous differential equation (see 
Ref. 1121 for more details). The corresponding homogeneous differential equation has 
two solutions J\{u) and J2(u). In order to reduce the transcendentality we may 
substitute Szlu) with J\(u) or J-i[u). The analytical expressions of J\ and J2 are 
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given in the Appendix of Ref. [T2j for example, if < u < 1 they read 

Ji(u) = 1 K(a(u)) , 

V'(l + ^) 3 (3-V^) y ' 

J 2 (u) = 1 K(l - a(u)) , 

V(i + V^) 3 (3-V^) 

(1- v^) 3 (3 + V^) 
(1 + V u ) ( 3 ~ V u ) 
In Ref. [12] we found also that 

-1 

du Jx(u) = Cl 2 (7r/3) . (59) 



We consider here the integrals of products of J; : we have checked up to 2400-digits 
of precision the following equalities: 



f 1 1 

/ du Jf(u) = -A 4 , 

Jo 8 



(60) 



d« 4{u) = -A 4 , (61) 



f 1 1 

/ d-u Ji(m)J 2 (m) = -£? 4 , 

Jo 4 



(62) 



and 

r9 



3 

Jl{u) = -A4 , (63) 



9 



Ja(w) = -^4 , (64) 



,0 j 
J du J 1 (u)J 2 (u) = -B 4 



(65) 



Therefore we have identified some one-dimensional integral representations of the 
numerical constants A 4 and B 4 . 

3.5. The key observation 



Eqs. (j60|) - (j65j) are not satisfactory elementary definitions of A 4 and B 4 , because of 
the complexity of the arguments of the K function in Eg. ([58]) . Our aim is to find 
out integral representations of A 4 and B 4 as simple as Eq. (|42"|) and Eq. ([45|) . We 
make a comparison between the integrals of products of K of the family (|40|) and 
the integrals of products of Ji. Eq. ([59|) corresponds to 



1 , K(t) ir 2 
dt , v ' = — 66 

We try to modify the integrand of so that the result contains the constant 
C12(tt/3). Our many year experience with the analytical calculation of 3-loop g-2 
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suggests that such a constant is usually associated with integrals containing the 
polynomial 1 + 1 + t 2 in the denominator. This is a factor of 1 — t 3 . Therefore we try 
to consider a "cubic" modification of the usual elliptic integral K(m). One fruitful 
choice is 

KJm) = / ; , EJm) = / — J, , (67) 

V ' Jo ^/(l-t 3 )(l-mt 3 ) 2 V ' Jo VT^¥ V ' 

or, equivalently, expressing them in terms of the hypergeometric function 



# c (m) = -|= 2 il(H;m) , E c (m) = ^ 2 F 1 ^-i;m) . (68) 



Now we calculate the numerical values of the integrals obtained by replacing K with 
K c in Eq. (|4ll)) . and we look for relations with the constants A 4 and B 4 . Luckily, we 
find 

M = \[ X d* K M™L- X \ (69) 
B^[\ X ^M. (70) 

4 J v ' 

We stress the tremendous simplification obtained by going from the usual descrip- 
tion with elliptic integrals l(gg |) -l|g5 ]) to the "cubic" version (f68|l - (f70j) . 



4. Four-loop results 

For the sake of brevity we define the following constants 
rl K 2 (x) 

C = I dx ; cW = 7.396 099 534 768 919 553 449 114 417 961 526 519 642. 
V 1 - x 



v= r d x K ^ K ^- x) 



(71) 



4.860 756 383 778 595 063 430 474 772 965 586 029 529. . . , 

(72) 

1-1 E 2 (x) 

E = I dx ; cV ' = 2.376 887 326 184 666 003 152 855 958 761 330 926 023 ... . 

(73) 



o 



Now we look for relations between the numerical values of the above constants C, T> 
and E and the numerical values of S 4 (D = 2) Eq.(|28|l and S±(D = 4 - 2e) Eq.([29]). 
We find that 

S A {-1, 1,1,1,1,1,1)= 2) =7rV3C , (74) 

or, alternatively, 

S 4 (-l, 1,1, 1, !,!,£>= 2) = ^nB 4 ■ (75) 
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we note the appearance of the factor 4n/3, similar to the appearance of 47r/-\/l5 in 
Eq.([7]). By using the integer-relation search program PSLQ^we have been able to 
fit the numerical result of Eq. ([2"9")l with the analytical expression 

5 45 4255 106147 

ttV3 , , 2320981 

C ^W( 297C - 1477 ^-^0736-- 

(76) 



The equalities Eg. ([74)1 and Eq. (|76| are the main result of this paper; they have 
been checked up to 2400 digits of precision. Note that constant V does not appear 
in Eq.jH and Eq.j76|). 



5. Four-dimensional lattice integrals 

Considering lattice perturbation theory, at one loop level one finds these integral!^ 
f n d 4 k 1 

Za= 7— — rr 5 = 0.154 933 390 231 060 214..., (77) 

J_ K (2^) 4 4EL 1 sin 2 (fc A /2) 1 ' 

and 

= r ^ 4 fc ^ 2 (fci/2)sin 2 (fc 2 /2) = ^ 

Some years ago, while visiting the Department of Physics of Parma, York Schroder 
pointed out to the author that, whether 3-loop g-2 was known in analytical form, 
no analytical result was known for the "simple" lattice 1-loop tadpole Zq. Puzzled 
by this fact, and noting that Zq can be reduced to a triple elliptic integral, we have 
tried to relate the numerical values of Zq and Z\ to the new constants C, T> and 
E. Working with only 10-digits precision numbers we have discovered numerically 
that 



S 4 (-l, 1,1, 1,1, 1,^=2) 



That is 



8/3 . (79) 



Zott" = . (80) 

By using again PSLQ, we have also found that 

Z^ = -^(ZC + 7E) + ^-l. (81) 

Values of Zq and Z\ with 400 digits of precision are quoted in Ref. [13 Very kindly, 
York Schroder provided us 16000-digits values. By using these numbers, we have 
checked Eq. ([80]) and Eq. ([8Tj) up to 2400 digits of precision, the maximum precision 
of our values of C and E. 
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We note also that the integral ([77)) can be rewritten into the so-called Watson 
integral in 4-dimensions (see Refs. [16j ITT)) 

uU) = — f f f r — ^ dh ^ = 8 Z Q (82) 

(2tt) 4 J_ n J_ n J_ 7T J_ 7T 4 - cos ki - cos k 2 - cos fc 3 - cos fc 4 ' 

From Eq. ([8T)]) one obtains 

m(4)tt 3 = 3V3C . (83) 



6. Conclusions 

In this paper we have at last identified beyond any reasonable doubt the analytical 
constants which appear in the simplest non-trivial 4-loop g-2 master integral. We 
have also discovered that the same constants appear in some 4-dimensional lattice 
integrals. Clearly we still do not know a rigorous proof of these relations, but, once 
the form of the re sults i s known, we hope that proofs will be easier to find (see the 
very recent paper J33H1). 



6.1. Note on Ref. \M 

While we were completing this paper, kindly David Broadhurst sent us a copy of his 
new papeJ^. In that paper, several elliptic integral evaluations of Bessel moments 
are performed. In particular, a proof of our Eq.© and Eq. ([22|) is given, as well as 
of Eq. ([44|) and Eq. ([46|) . Constants analogous to our A and B are found, 04,0 = 27ryl 
and S4.0 = B, as well as ^6,1 — A4/8 and Se,i — S^{D = 2)/16. In addition, several 
relations between Bessel moments are found, and some evaluations of double elliptic 
integrals are done. 
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